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I.  Introduction 


This  report  summarizes  the  research  performed  under  DAAG-29- 

81-K-0175.  The  objectives  of  the  research  were  to  assess  the 

merits  of  resonant  frequency  doubling  in  paramagnetic  spin  systems 

in  the  near  millimeter  (NMM)  region.  The  motivation  for  this 

occurs  because  the  efficiency  of  conventional  solid  state  devices 
_  2 

falls  as  oi  so  that  a  doubler  with  a  conversion  efficiency  of 
>  25%  would  have  a  higher  system  efficiency  at  2w  than  a  primed 
power  source. 

In  Sections  II  and  III  we  review  the  main  characteristics  of 
resonant  frequency  doubling  from  an  optical  standpoint  in  a  three- 
level  system.  In  Section  IV  we  consider  the  characteristics  of 
available  spin  ions  in  candidate  host  materials  and  show  that 
the  approach  does  not  seem  promising  because  of  inordinately  long 
crystal  lengths.  In  Section  V  we  comment  on  resonant  frequency 
tripling  in  a  three- level  system  and  in  Section  VI  we  comment  on 
the  possibilities  which  might  exist  in  heterostructure  super¬ 
lattices. 


II.  Overview 


In  an  earlier  paper1  a  semiclassical  treatment  of  second  har¬ 
monic  generation  in  a  three-level  system  was  considered,  motivated 

by  recent  experimental  observations  of  doubling  in  iron  doped  rutile 

2 

in  the  microwave.  A  density  matrix  description  of  three-wave 
parametric  interactions  was  used  to  obtain  analytical  nonperturba- 
tive  solutions  for  the  nonlinear  polarization  within  the  limits  of 
the  rotating  wave  approximation,  RWA.  In  Fig.  1  the  perturbation 
results  are  presented  for  the  second  harmonic  polarization  in  a 
three- level  system  with  and  without  the  RWA.  Clearly,  this  is  an 
excellent  approximation  for  the  nearly  resonant  and  resonant  cases. 

The  nonlinear  polarizations  were  next  coupled  with  Maxwell's 
equations  and  solved  numerically.  For  a  strongly  saturating  pump 
and  for  certain  matrix  elements,  large  conversion  efficiencies 
(>  50%)  were  predicted.  These  results  were  ascribed  to  the  large 
nonlinear  polarization  associated  with  a  resonant  system  and  to 
AC  Stark  splitting  of  the  harmonic  transition  for  >>  A^2  which 

minimizes  absorption  of  the  harmonic.  Additionally  the  pump  decay 
is  linear  with  distance. 

In  this  report  we  first  expand  on  the  results  presented 
earlier  for  harmonic  generation  in  a  three-level  system  and  then 
we  apply  these  results  to  magnetic  resonance  systems.  Clogston3 
was  the  first  to  treat  this  problem  albeit  in  the  limit  of  a 
strong  pump  and  a  weak  harmonic.  Related  are  calculations  of 
frequency  conversion  in  a  four-level  system  again  in  the  strong 


Figure  1.  Perturbation  calculation  of  SHG  polarization  with  and  without  the  rotating 
wave  approximation,  ft, _  =  =  50  GHz. 
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pump  limit,  and  calculations  of  harmonic  generation  m  two  level 

5  6  7 

spin  systems  '  and  in  fully  nonresonant  systems. 

In  addition,  many  authors  have  reported  experimental  obser- 

Q 

vation  of  nonlinear  effects  in  paramagnetic  materials.  Kellington 

was  the  first  to  observe  second  harmonic  generation  from  ruby 

5  6  9 

crystals  at  room  temperature.  Boscaino,  et  al.  '  '  extended  this 

work  to  include  low  temperatures  and  saturation  effects.  A  number 

of  other  experiments  have  been  reported  in  resonant^"® '  ^  and 

either  nonresonant  or  two-photon  resonant  systems.  '  The  best 

-4 

conversion  efficiency  (10  )  was  obtained  by  Yngvesson  and 

2  +3 

Kollberg  in  Fe  /rutile  at  liquid  helium  temperature  for  the 
fully  resonant  case.  Encouraged  by  these  results,  we  calculate 
the  practical  limits  of  resonant  second  harmonic  in  this  system. 
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III.  General  Considerations 
a)  Basic  Equations 

In  this  section  some  definitions  and  results  presented 
in  an  earlier  paper^  and  relevant  to  the  following  are  reviewed. 

The  initial  growth  of  the  harmonic  intensity  accompanying 
the  linear  decay  of  a  saturating  pump  was  shown  to  be 

cR.r  , 
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sionless  space  variable  and  p^a  is  the  parametric  contribution  to 
the  second  harmonic  intensity.  The  latter,  again  in  the  strong 
pump  limit, can  be  written  as  a  function  of  the  saturated  population 
differences , which  may  be  positive,  negative  or  even  zero,  or 
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complex  detunings,  c,  the  space  variable,  is  related  to  an  actual 
propagation  distance  by  the  line  center  Beer's  coefficient,  <*2.2' 

ct2  - 


for  the  |1>  -*•  1 2>  transition. 


t 


a12  z/2  therefore,  if  the 
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Rabi  frequencies  are  expressed  in  units  of  T2,  C  measures  the 
distance  in  units  of  field  Beer's  coefficient.  The  harmonic 
intensity  is  normalized  to  the  pump  intensity  by  defining 
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b)  Temperature  Dependence 

A  particularly  simple  limit  for  p 
when  T^=T2  and  all  matrix  elements  are  equal. 


a 
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is  arrived  at 
Then 
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(4) 


which  implies  an  interference  effect  as  |2>  is  populated  and  hence 
a  decrease  in  the  value  of  p^  with  an  increase  in  temperature. 
Therefore,  from  Eq.  (1)  the  initial  growth  of  the  harmonic  inten¬ 
sity  requires  a  larger  conversion  distance.  No  significant  change 
in  the  peak  conversion  efficiency  is  predicted  in  a  three-level 
system,  however,  if  other  nearby  levels  are  populated  a  reduction 
in  the  harmonic  intensity  would  be  expected.  In  Fig.  3  the  peak 
conversion  efficiency  and  distance  are  plotted  for  T1=T2  and  Ti  = 
100  T»  for  A, _T_  =  4  and  R  =  6.  In  each  case  for  AE/kT  >>  1  the 

2  12  2  U 

conversion  distance  is  observed  to  increase  sharply  with  tempera¬ 
ture.  No  significant  change  in  the  peak  efficiency  is  predicted. 

c)  Tx  ±  T: 

The  effect  of  T^T2  can  be  qualitatively  predicted 

14 

from  the  T^  dependence  of  the  population  terms. 


CONVERSION  DISTANCE 
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These  equations  are  derived  by  considering  the  population  equa¬ 
tions  for  a  resonant  field  and  a  negligible  A.-^ 3 .  With  an  increase 
in  T^  the  population  differences  decrease  and  therefore  satura¬ 
tion  occurs  more  quickly,  the  pump  decay  is  slower  and  from  Eq. 

(1)  the  build  up  of  the  harmonic  is  slower  and  more  spread  out. 
Figures  4  and  5  are  plots  of  maximum  conversion  efficiency  and 
distance  for  several  T^/T2  rat^os  anc*  a  variety  of  matrix  ele¬ 
ment  ratios.  The  conversion  distance  is  predicted  to  increase 
sharply  while  little  change  in  peak  conversion  efficiency  is 
expected.  Previously,  for  T^=T2  a  resonance  was  observed  for 
R  =  1//2.  A  similar  resonance  is  not  calculated  when  T-^/T^ 

d)  Variation  of  Individual  Transition  Probabilities 

Transition  probabilities  between  different  pairs  can 
differ  significantly  depending  on  the  forbiddance  of  the  transi¬ 
tion  between  states  and  on  the  degree  of  mixing  of  states.  The 
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Figure  5 


Conversion  efficiencies  and  conversion  distances  (in 
relative  units)  as  a  function  of  ^2T2  f°r  x*y=*z*0 

25  (increments  of  .25) 


R,  *.25-1. 

u 


*12T2 
for  T^=l *10 


CONVERSION  EFFICIENCY  I^x100(%) 


CONVERSION  DISTANCE 
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general  form  of  the  population  equations  can  be  written' 


fH'.Pljj  +  l 
k 


(pkka,kj_0iiwjk) 


(6) 


where  is  a  per  unit  time  transition  probability  subject  to 

the  constraints  of  detailed  balancing.  A  relaxation  time  can  be 

i 

defined  as  TV ^  =  T\  ^  P^j.  In  Figs.  6  and  7  we  show  the 

effect  on  the  conversion  distance  of  decreasing  respectively  T^2 
and  T13  from  100xT2  to  T2.  All  other  T  —  are  set  equal  to  100* 
T2-  In  each  case  as  T^2  or  T^3  approach  T2  the  conversion  dis¬ 
tance  approaches  that  calculated  for  a  decrease  in  the  overall 
T^.  However,  a  decrease  in  T23  does  not  significantly  affect 
the  conversion  distance  as  is  expected  considering  the  interfer¬ 
ence  effect  associated  with  populating  \2>. 


e)  Effect  of  Detuning 

In  Figs.  8,  9,  and  10  the  conversion  efficiency  and 
conversion  distance  are  plotted  versus  A12T2  for  a  system  where 
the  pump  is  detuned  by  10  linewidths  and  the  harmonic  is  resonant 
and  for  the  case  where  the  pump  is  detuned  by  10  linewidths  and 
the  harmonic  is  detuned  by  20  linewidths.  In  each  case  a  decrease 
in  the  conversion  efficiency  is  predicted  as  well  as  in  increase 
in  the  conversion  distance.  Both  the  conversion  efficiency  and 
conversion  distance  converge  to  the  resonant  case  as  the  line  is 
power-broadened . 

These  results  are  consistent  with  the  non-zero  absorption  which 
is  observed  when  the  system  is  not  on  one  and  two-photon  resonance.* 


CONVERSION  EFFICIENCY  I^x100(%) 
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Figure  7 


Variation  of  conversion  efficiencies  and  conversion 
distances  (in  relative  units)  with  T12  *  T23  * 

100><T2,  x»y»z*0,  *  1  and  *  1-6. 
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Figure  9.  Plot  of  conversion  efficiencies  and  conversion  dis¬ 
tances  (in  relative  units  for  x=-y*10  and  z=0,  p?,  *1 
and  Ry  =  l-6.  T^lOOxT^. 


CONVERSION  DISTANCE 
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TAz 


Figure  10.  Plot  of  conversion  efficinecies  and  conversion  dis¬ 
tances  (in  relative  units)  for  x=y*10,  z*20,  p?.  *1 

and  R  = 1-6. 

M 
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f )  Effect  of  Inhomogeneous  Line  Broadening 

A  treatment  of  inhomogeneous  line  broadening  follows 
directly  from  the  consideration  of  detuning  in  the  previous  sec¬ 
tion.  If  we  consider  a  system  with  a  Gaussian  distribution  of 
detunings,  for  example,  the  velocity  distribution  in  a  Doppler 
broadened  line,  then  the  total  nonlinear  polarization  is  the  sum 
of  the  nonlinear  polarization  associated  with  each  velocity  point 
(detuning) .  In  Figs.  11  a)  and  b)  peak  conversion  efficiencies 
and  conversion  distances  are  plotted  for  several  T^/T2  ratios  and 
a  variety  of  inhomogeneous  linewidths.  As  expected  for  inhomo¬ 
geneous  linewidths  which  are  larger  than  T2  the  loss  in  conver¬ 
sion  efficiency  is  significant  (>  25%)  as  is  the  increase  in  con¬ 
version  distance.  Again  the  results  converge  to  the  homogeneous 
case  as  the  line  is  power  broadened.  In  general,  as  for  the  homo¬ 
geneous  case,  the  larger  matrix  elements  result  in  the  best  con¬ 


version  efficiencies. 
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IV.  Application  to  Magnetic  Resonance  Systems 

2 

Encouraged  by  the  report  of  Yngvesson  and  Kollberg  of  ef¬ 
ficient  second  harmonic  generation  in  iron  doped  rutile  we  have 
calculated  the  practical  limits  of  conversion  efficiency  in  this 
and  similar  paramagnetic  materials.  Materials  which  exhibit  nuc¬ 
lear  magnetic  resonance  are  not  considered  since  the  moments  are 
a  thousand  times  less  and  the  energy  level  splitting  are  much 
smaller. 15 

Iron  doped  rutile  is  a  convenient  system  to  use  as  a  model 
since  much  is  known  about  its  magnetic  resonance  and  material 
properties  due  to  an  earlier  interest  in  maser  materials. 17-1^ 

In  addition  it  is  representative  of  a  large  group  of  paramagnetic 
compounds  the  transition  metals  and  the  loss  tangent  is  relatively 
small  (£  10  5) .  In  rutile,  the  Fe+3  ion  is  situated  at  two  magnet¬ 
ically  inequivalent  sites  of  Djh  symmetry  which  are  related  to 
each  other  and  to  the  optic  axis  by  a  rotation  of  90°.  The  spin 

Hamiltonian  is  characteristic  of  an  ion  (s=5/2)  at  a  site  of 

20  21 

tetragonal  symmetry.  '  The  spin  Hamiltonian  parameters  have 

22 

been  shown  to  be: 

D  =  20.35  GHz  E  =  2. 21±0.07  GHz 

a  *  1.1  to  0.6  GHz  F  =  -0.5±0.3  GHz 

Inasmuch  as  the  zero  field  splitting  is  large,  substantial  energy 

level  splittings  are  obtainable  with  small  magnetic  fields.  Spectra 

of  the  two  inequivalent  sites  coincide  in  a  plane  defined  by  the 
c-axis.  Figure  2  is  ex  energy  level  diagram  for  Fe  in  rutile  at 


21 


an  angle  45°  from  an  a-axis  and  32. 5° from  the  a-a  plane.  A  possi¬ 
ble  parametric  interaction  is  indicated.  In  this  system  at  a  con- 

20  -3  —8 

centration  of  1*10  cm  and  an  estimated  T 2  of  lx10  sec,  the 

scale  factor,  a12/2'  aPProx;i-ma'tely  15  cm-1  at  a  frequency  of 

3  cm-’'".  It  should  be  noted  that  concentrations  of  this  order  of 

23 

magnitude  are  unrealistic  in  rutile  although  realizable  in  hosts 
(spinels,  corundum,  etc.)  where  problems  with  charge  compensation 
do  not  occur. 

The  longitudinal,  ,  relaxation  properties  of  Fe+^  doped  rutile 

24—28  28  20 

are  well  characterized.  ’  '  Above  10°K  T^  is  dominated  by 

Raman-type  processes  and  at  room  temperatures  the  relaxation  time 

is  controlled  by  direct  processes  and  is  on  the  order  of  a  few 
25 

milliseconds.  This  is  roughly  four  to  five  orders  of  magnitude 
larger  than  T2«  Neglecting  for  the  moment  inhomogeneous  line 
broadening  and  assuming  only  the  ground  state  to  be  populated, 

Fig.  4  d)  is  representative  of  Fe+^/rutile  at  liquid  helium  tem¬ 
peratures.  From  this  plot  and  considering  a  scale  factor  of  15 
cm  ^  we  see  that  for  conversion  efficiencies  on  the  order  of  25% 
or  better  the  crystal  would  need  to  be  roughly  a  thousand  centi¬ 
meters  in  length.  From  Fig.  3  we  see  that  a  further  increase  in 
conversion  distance  is  expected  if  we  consider  that  AE/kT  ^  1 
and  therefore  approximately  one-third  of  the  total  population  is 
not  in  the  ground  state. 

Apparently,  in  order  to  achieve  realistic  conversion  distances, 
T^  must  be  decreased  relative  to  T2.  At  room  tempeature  T^ 


approaches  but  all  levels  become  almost  equally  populated. 
Referring  again  to  Fig.  3  for  AE/kT  =  .01  and  T1=T2  a  hun<*re<* 
centimeter  crystal  would  be  required.  Another  possibility  is  to 
exploit  nonresonant  and  resonant  cross-relaxation  processes  between 
Fe+3  and  faster  relaxing  ions.  Cross-relaxation  has  been  exten¬ 
sively  investigated  in  ruby  and  is  attributed  to  magnetic  dipole 

and  virtual  phonon  interactions  as  well  as  quadrapole  electric 
31-33 

field  effects.  Of  course,  the  maximum  decrease  in  is 

observed  for  resonant  processes  where  an  impurity  would  contribute 

pump  absorption.  can  also  be  decreased  by  increasing  the  con- 

27 

centratron  of  ions  in  the  host.  This  is  probably  related  to 

an  increase  in  exchange  interactions.  The  conversion  distance 

can  also  be  decreased  by  increasing  T2.  It  has  been  shown34 ' 5 ' 6 ' 9 ' 35 

that  when  the  oscillating  magnetic  field  is  much  larger  than  the 

local  field  seen  by  each  ion,  T2  can  be  effectively  reduced  for 

the  pump  transition.  This  has  been  demonstrated  in  ruby  and 
5  6  36 

DPPH,  '  '  and  more  recently  for  an  optical  transition  in 
5  2 

Pr :LaF 3 .  In  the  former  case,  TL  was  reduced  to  10~7  sec.  It 

is  limited  by  local  fields  with  a  finite  correlation  time.3^  in 

ruby  the  upper  limit  on  T2  is  determined  by  surrounding  ^7A1 

ions  (1=5/2)  of  the  host  lattice.  In  rutile  roughly  90%  of  the 

titanium  sites  are  occupied  by  isotopes  with  1=0.  Therefore,  a 

somewhat  longer  T2  might  be  realized.  In  addition  an  anomolously 

low  T2  has  been  observed  for  transitions  close  to  the  zero  field 
37  38 

splitting.  '  It  is  also  important  to  remember  that  increasing 
T2  also  increases  the  scale  factor  hence  further  reducing  the 
actual  material  length. 
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If,  in  fact  can  be  made  roughly  equal  to  T2  and  T2  can  be 
reduced  to  10~®-10-7  sec  substantial  conversion  efficiencies  at 
moderate  conversion  distance  might  be  feasible.  As  an  example, 
for  T2=T^=1x10-^  sec  the  scale  factor,  assuming  a  concentration 
of  1*10^  cm-"*  is  150  cm-^.  For  A^2T2  =  10  the  conversion  distance 
is  on  the  order  of  0.2  cm  and  for  certain  combinations  of  matrix 
elements  will  approach  100%. 

Of  course  random  crystal  defects  and  impurities  will  result  in 

39 

local  field  variations  from  site  to  site  in  the  host.  Therefore, 
the  linewidths  of  magnetic  resonance  systems  are  usually  inhomo- 

*  _7  _  q 

geneously  broadened  with  T2  varying  from  10  -10  sec.  This 

effect  enters  into  the  calculation  in  two  ways.  First,  the  scale 

* 

factor  is  now  calculated  from  T2>  the  inhomogeneous  linewidth. 
Second,  a  spin  system  must  now  be  considered  where  each  individual 
spin  packet  is  not  on  resonance  but  rather  a  Gaussian  distribution 
about  full  resonance  is  assumed. 

Linewidths  observed  for  iron-doped  rutile  at  high  concentrations 
19  -3 

(1x10  cm  )  are  very  broad,  several  hundred  Gauss  at  x-band 

_  .  23  *  _q 

frequencies  corresponding  to  a  T2  of  <1x10  sec  and  hence  a 

scale  factor  of  <.15  cm”^.  With  Tj,=T2  an<^  only  the  ground  level 

populated  a  crystal  30  cm  long  would  still  be  required. 

Unusually  narrow  linewidths  as  well  as  short  relaxation  times 

generally  characterize  lanthanide  ion  spectra  but  energy  level 

-1  40 

splittings  usually  range  from  10-100  cm  .  Some  notable  excep¬ 
tions  are  Dy+^  and  Er+^  in  nearly  cubic  hosts  like  YAgO^  where 


2 


4 

energy  level  splittings  can  be  on  the  order  of  1  cm.  Unfortun¬ 
ately,  observed  linewidths  are  very  broad. 

Non-Kramers  ions  for  example,  Fe  (T^=2  us)  and  ions  which 
exhibit  Jahn-Teller  coupling  to  the  lattice,  V+^,  (T^=l*10  ^  sec)^ 

might  be  considered  but  again  the  inhomogeneous  linewidths  are  very 
wide,  in  general,  for  good  conversion  efficiencies  at  moderate 

conversion  distances  it  is  necessary  that  be  of  the  same  order 

* 

of  magnitude  as  T2  and  neither  T2  nor  T2  should  be  faster  than 
10  or  10  sec.  In  addition,  the  bulk  of  the  equilibrium  popu¬ 
lation  must  be  in  the  ground  state  implying  for  NMM  wave  fre¬ 
quencies  liquid  helium  temperatures.  Some  other  candidate  magne- 

+  2 

tic  resonance  systems  might  be  Mn  in  CaW04  co-doped  with  a  fast 
+  2  +3 

relaxer  like  Fe  or  Ho  doped  YLiF^. 
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V.  Resonant  Tripling 

Because  of  the  material  parameter  limitations,  the  param¬ 
agnetic  spins  systems  do  not  appear  favorable  for  efficient  (>25%) 
second  harmonic  generation.  However,  the  tunability  and  choice 
of  energy  spacing  does  not  preclude  higher  order  harmonic  situa¬ 
tions  such  as  resonant  frequency  tripling. 

We  have  studied  this  problem  also  in  a  manner  analogous  to 
the  resonant  doubling  system  just  outlined.  The  equations  and 
resulting  solutions  are  totally  opaque  and  will  not  be  listed 
for  this  reason. 

Our  findings  on  the  dependence  of  the  third  harmonic  conver¬ 
sion  efficiency  versus  various  material  parameters  and  pump  in¬ 
tensity  suggest  that  unlike  the  case  of  second  harmonic  generation 
in  a  three- level  system  for  which  the  harmonic  absorption  is  split 
out  of  resonance  by  the  pump,  harmonic  absorption  is  split  out  of 
resonance  by  the  pump,  harmonic  absorption  is  always  present  and 
limits  the  conversion  efficiency  to  about  ^  25%.  Further,  under 
conditions  of  strong  conversion  into  the  third  harmonic,  there 
is  a  non-negligible  gain  for  the  second  harmonic  which  places 
some  severe  restraints  on  the  phase  matching  system  to  act  as  a 
filter  for  this  work.  Because  of  these  features,  resonant  tripling 
in  general  does  not  appear  to  be  a  viable  approach  if  high  effici¬ 
ency  is  sought. 
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VI.  Nonlinear  Semiconductors 

In  an  effort  to  find  materials  more  suitable  for  efficient 
doubling,  some  thought  was  given  to  semiconductor  systems.  It  is 
well  known  that  the  I-V  curves  of  a  diode  are  highly  nonlinear 
with  the  forward  biased  situation  used  in  mixer  technology  and 
the  reverse  biased  situation  used  in  varactor  technology.  In 
the  latter  case,  the  nonlinearity  arises  because  of  the  dependence 
of  the  depletion  width  and  hence  depletion  layer  dipole  on  the 
applied  potential.  The  nonlinearity  is  enormous  by  optical  stan¬ 
dards  but  is  more  distributed  and  requires  conduction  paths  to 
transform  free  space  or  guided  waves  into  the  devices.  Out  of 

necessity,  this  introduces  losses  so  that  the  best  reported  doubl- 

44 

ing  efficiency  in  the  NMMW  was  'v  10%.  Other  schemes  involving 
modulation  doped  heterostructures  might  be  anticipated  to  be  no 

4  5 

better  although  still  worthy  of  investigation. 

The  primary  origin  of  the  optical  nonlinearity  in  bulk  semi¬ 
conductors  is  the  non-parabolicity  of  the  bands.  While  of  modest 
size  in  the  infrared,  the  nonlinear  coefficients  become  very 
large  a  NMMW  frequencies.  Even  with  this  increase,  the  required 
power  for  efficient  conversion  is  still  large,  'v  MW/cm  for  a  corn- 
version  ^  50%  in  1  cm. 

Highly  nonparabolic  bands  are  known  to  exist  in  heterostructure 
systems  comprised  of,  for  example,  GaAs-A^Ga^^As  superlattices. 
Tsu  and  Esaki^  were  the  first  to  consider  optical  nonlinearities 
,  in  a  symmetric  periodic  super lattice.  They  showed  that  with  a 
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symmetry  breaking  DC  field,  the  AC  current  has  a  component  at 

2oj  corresponding  to  doubling  wherein  the  nonlinearity  is  due  to 

47 

an  I-V  nonlinearity.  More  recently  Bloss  and  Friedman  and 
48 

Yuen  have  both  considered  third  order  optical  nonlinearities 
in  superlattices  and  showed  that  these  may  be  large  under  suit¬ 
able  conditions,  some  100  times  larger  than  in  the  bulk. 

The  above  three  theoretical  studies  were  all  based  on  non¬ 
linear  transport  properties  of  carriers,  located  in  the  lowest  sub¬ 
band  within  the  conduction  bands.  It  is  also  possible  to  have 
radiative  transitions  between  subbands  and  to  have  optical  non- 
linearities  in  analogy  with  what  is  seen  in  normal  situations. 

Orlov  has  considered  a  symmetric  superlattice  and  concluded  that 

the  third  order  optical-like  nonlinearity  may  be  some  six  orders 

49 

of  magnitude  larger  than  in  the  equivalent  bulk  crystal.  The 
origin  of  the  nonlinearity  is  in  the  large  optical  matrix  elements 
associated  with  the  envelope  wave  functions  of  the  confined  car¬ 
riers.  This  is  analogous  to  the  situation  in  the  nonlinear 
polymer  MNA. 50 

For  symmetry  reasons,  the  lowest  order  nonlinearity  associated 
with  a  periodic  symmetric  potential  is  restricted  to  odd  har¬ 
monic.  We  have  recently  shown  that  asymmetric  potentials  may  be 
achieved  in  superlattices  with  the  use  of  compositional  grading 
and  that  the  lowest  order  quadratic  nonlinearity  is  estimated  to 
be  some  two  orders  of  magnitude  larger  than  in  bulk  samples. ^ 
Because  of  practical  difficulties  in  growing  very  low  compositional 
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gradients,  this  approach  will  be  most  useful  in  the  infrared 
region  and  is  expected  to  yield  nonlinearities  in  the  NMM  region 
comparable  to  existing  materials. 

There  are  other  possiblities  for  the  superlattice  case  in¬ 
volving  charge  rearrangement  within  a  potential  region,  analogous 

45 

to  the  varactor  case,  which  have  not  yet  been  explained  in  any 
detail.  These  are  not  radiative  transitions  based  nor  nonlinear 
I-V  in  nature  but  rather  "loose"  dipoles.  In  the  NMM  region  some 
power  dissipation  must  occur  out  of  necessity  but  the  magnitude 
of  the  loss  and  the  magnitude  of  the  nonlinearity  remain  to  be 
determined.  Clearly,  there  are  other  possibilities  for  efficient 
doubling  in  the  NMM  range  and  some  of  these  approaches  seem 
worthy  of  further  study. 
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